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We analytically examine the dynamics of quasi-particle 
modes occuring in a Bose-Einstein condensate which is sub- 
ject to a weak acceleration. It is shown that the momentum 
of a quasi-particle mode is squeezed rather than accelerated. 



The achievement of Bose-Einstein condensation in 
magnetic traps |j| has led to a great interest in the phys- 
ical properties of inhomogeneous atomic quantum gases. 
In the current research, one of the central topics is to 
study collective excitations of quasi-particle modes in 
Bose-Einstein condensates (BEC). The frequency spec- 
trum of elementary excitations in a trapped 87 Rb BEC 
has been experimentally measured at JILA [^j. The 
propagation of sound waves in an excited BEC has been 
demonstrated at MIT j|. Theoretical studies have re- 
vealed the characteristics of low-energy quasi-particle 
modes in such nonlinear systems Q. 

In this paper we are interested in another aspect of 
collective excitations of BECs. We will study the dy- 
namics of quasi-particle modes and show theoretically 
how they are manipulated by a weak homogeneous force. 
Since the dispersion relation of quasi-particle modes is 
much different from that of ordinary Schrodinger parti- 
cles their dynamics will show up unusual effects like an 
acceleration-induced squeezing of position or momentum. 
A squeezing of the position may correspond to a local- 
ization of sound waves in quantum gases. 

To initiate our idea, we start with a simplified model 
where a BEC trapped in a box of size L. In practice the 
box may be realized by a sufficiently anharmonic trap or 
an atomic cavity. We then add to this trapping potential 
a second, linear potential which describes the external ho- 
mogeneous force. As a result, the complete potential be- 
comes asymmetric and should incorporate an area where 
the potential is almost linear (see Fig. 1). A harmonic 
trap is not suited for this purpose since the homogeneous 
force then only leads to an overall displacement of the 
trap without changing its shape. 

The time evolution of a BEC can be appropriately de- 
scribed by the time-dependent Gross-Pitaevskii equation, 
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where ip(x) is the normalized collective wavefunction of 
the BEC, M the atomic mass, and k := infr 2 N a sca ± / M . 
N is the number of condensed atoms and a sca t > the 



scattering length of the interaction between the atoms. 
V(x) is an external potential. 

The ground state ip g (x,t) = exp[— ifit/h]tha(x) of the 
trapped BEC is a stationary solution of Eq. (|l|) , where [i 
denotes the chemical potential. ip g (x) can be assumed to 
be real. For current experiments the nonlinear interac- 
tion between the atoms is relatively strong compared to 
the BECs kinetic energy. It is then possible to neglect 
the kinetic term in Eq. (Q) and to describe the ground 
state in the Thomas-Fermi approximation (TFA) B , 
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Of course this approximation is only valid in areas where 
the r.h.s. (right-hand-side) does not become small or neg- 
ative. 

Elementary excitations around the ground state tp g of 
Eq. ([!]) are described by writing the collective wavefunc- 
tion as tj){x,t) = exp[— ifit/h]{t/> g (x) + Q(x,t)} and ex- 
panding Eq. (Q) up to terms linear in 0. The resulting 
equation, 
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with R := p 2 /(2M) + V - fx + 2n\^ g {x)\ 2 and S := Kip 2 , 
has for time-independent R and S the formal solution 
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(Go := <d(x,0)). The formal frequency operators 
h 2 oj± := y/(RTS)(R±S) 
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correspond to the spectrum elementary excitations which 
is usually calculated by using the Bogoliubov-DeGenncs 
equations (see, e.g., Ref. 0). We remark that one can 
rewrite Eq. (^) using or a>_ alone by exploiting rela- 
tions like (R — S) cos[ij^t] — cos[il> + t](R — S). 

In areas where the TFA is valid the operators uj± can 
be simplified to 
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and a similar expression for uj_ . If we restrict our consid- 
erations to quasi-particle modes with small momentum 
we can neglect the second kinetic term in this equation 
and arrive at 



1 



hui4 



y /p2 (c*-V(x)/M) 



(J) 



where we have defined the velocity of sound c := y/p,/M. 
Obviously this approximation corresponds to the restric- 
tion to the linear part of the dispersion relation, TiCj+ = 
c\p\, for elementary excitations in a homogeneous, free 
BEC (V(x) = 0). 

Though Eq. (Q) looks simple the square root on the 
r.h.s. is operator- valued and in general difficult to an- 
alyze. We therefore restrict our analysis to one spa- 
tial dimension and consider an accelerating potential 
V(x) = —Max with ah <C c 2 , where a is the acceler- 
ation and L the length of the condensate. We then can 
expand uj + in a, leading to 
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In the derivation we have used the commutation relation 
[x, sgn(p)] = 2ihVo, where Po is the projection operator 
on the state p = 0. 

Some caution is needed when expanding the square 
root of an operator as in Eq. (||). Already Eq. (Q) 
is only meaningful if c 2 — V/M — c 2 + ax is larger 
than zero, in agreement with the assumptions of the 
TFA. Hence our approach is only useful for wavepack- 
ets that are sufficiently localized inside the interval x € 
[— c 2 /a, — c 2 / a + L], In addition we have to demand that 
the quasi-particle momentum is not too small so that 
its De-Broglie- wavelength is not larger than L. This be- 
comes evident in the expansion (^) which diverges for 
p — > 0. Mathematically the inverse of an operator with 
eigenvalue zero cannot be defined. However, the inverse 
makes sense on a subspace that excludes this eigenvalue, 
and our restriction to sufficiently large p ensures that we 
are working in this subspace. Our approach therefore al- 
lows us to study wavepackets which are localized in phase 
space. For simplicity we will assume that the wavepacket 
is localized around some momentum pq > in momen- 
tum space so that we can set the operator sgn(p) equal to 
one. This causes no loss of generality as any wavepacket 
subject to the constraints given above can be linearly 
split into a part with p > and one with p < 0. The 
latter case can be treated analogously to the case p > 
which we will consider now. 

Having derived an explicit expression for the frequency 
operators we can return to Eq. (^) and consider the evolu- 
tion of a localized elementary excitation. Obviously the 
time dependence of O is essentially determined by the 
four exponentials exp[±icD±t]. Since iD_ is very similar 



to 0)4- we will only study the exponentials of lu + which 
can be decomposed as 
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where the exponentials are correct to order a 2 . 
The proof of this relation is analogous to that 
of the Baker-Campbell-Hausdorf theorem. For in- 
stance, an expression like F(t) := exp[it(Xp + X'px 2 )] 
can be decomposed by making the ansatz F(t) — 
exp[iu(t)p] exp[iv(t)px] e?q>[iw(t)px 2 ] and differentiating 
both expressions for F(t) with respect to t. Using sim- 
ple commutation relations the resulting equation corre- 
sponds to a set of ordinary differential equations for u(t), 
v(t), and w(t) which can be solved exactly or, as is suffi- 
cient for our purposes, to order o 2 . 

Eq. (|9|) allows us to study in detail the time evolution 
of a quasi-particle mode Q(x,t) by discussing the action 
of exp[±iw + t] on |0 O ) 0, with {x\Q Q ) := 8o(ir). Before 
doing so we remark that for a = we have u)_ = cj+ 
and exp[±iLj + t] = exp[±ipct / h] (in accordance with the 
discussion given above p > has been assumed). The 
r.h.s. of this expression corresponds to a shift of the 
position, x' = x ± ct. Hence, for a = the exponentials 
exp[±iw + i] describe the propagation of sound waves to 
the right or to the left, respectively. As the right-moving 
part can be obtained from the left-moving part by time- 
inversion we will discuss the -left moving part ®l (x,t) :— 
(x\ exp[— wI>_|_i]|©o) only. 




FIG. 1. Sketch of the evolution of a localized quasi-particle 
mode in the linear regime of the combined potential V{x). 
A small excitation Oo of the ground state ip g is splitted into 
two counterpropagating sound waves. The homogeneous force 
induces a squeezing operation on these two wavepackets. The 
additional potential provided by ipg leads to a reduction of 
the acceleration induced by the force. 

To get an explicit expression for ©l(x, t) we have to an- 
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alyze the action of each individual exponential of Eq. (||). 
The first one is just a c-number and corresponds to a 
rescaling of the wavefunction (note that Eq. (||) does not 
conserve the norm of 0). As discussed above the sec- 
ond exponential corresponds to a time-dependent shift 
of the wavefunction and describes an accelerated sound 
propagation: 
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In this equation \<f>) denotes an arbitrary element of the 
Hilbert-space L 2 (R) of square integrable functions. The 
third exponential is of the form exp[— iapx] with a := 
at[l — at / (4c)] / (2hc) . This operator corresponds to a 
squeezing transformation, 
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A quick way to see this is to write the position and mo- 
mentum operator in terms of creation and annihilation 
operators as for the harmonic oscillator: x ~ a + a T , 
p ~ i(a' — a). The exponential then can be written as 
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1]} which is just the usual form of a 



single-mode squeezing transformation (see, e.g., Ref. 
Mathematically, Eq. ( |Tlj ) can be proven by differentiat- 
ing f(a,x) :— (x\e~ lapx \(f)) with respect to a. the re- 
sulting differential equation, df /da = —hd(xf)/dx has, 
for the inital value f(0,x) = (x\<fi), the solution (11). 
Analogously one can show that in momentum space the 
squeezing operator acts like 
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As usual the exponent of the squeezing factor exp[?La] has 
the opposite sign compared to the squeezing of the con- 
jugate variable x. New compared to the usual squeezing 
of ordinary particles is the c-number factor of exp[— ha] 
appearing on the r.h.s. of Eq. (|TT|) . This factor and 
the resulting asymmetry between Eqs. ([ll]) and (|lj) is 
a consequence of the operator ordering: px instead of 
(px + xp)/2. The order of px itself is a consequence of 
the asymmetric definition of lu + = [p 2 (c 2 + ax)] 1 / 2 . 

In an analogous treatment one can show that the 
fourth exponential, which is of the form exp[i/3pi 2 ] with 
(3 := a 2 t/(8hc 3 ), corresponds to a rational transforma- 
tion of x, 
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The singularity at x — l/(h/3) is unphysical since it only 
appears if the wavefunction is not localized inside the 
area where c 2 3> ax holds. The corresponding equation 
in momentum space is much more complicated and dif- 
ficult to interpret. We can gain some insight into this 



case by considering the Heisenberg-picture evolution of 
the momentum operator: 
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Since the momentum operator gets multiplied with the 
conjugate operator x its action in momentum space is 
difficult to interpret. However, it still remains an op- 
erator linear in p so that its action in position space 
is easy to understand: the wavefunction gets multiplied 
with (1 + h(3x) 2 and is then differentiated. 

The fifth exponential in Eq. (j|) acts as a position- 
dependent rescaling of the wavefunction. In mo- 
mentum space it corresponds to a imaginary shift 
of the momentum. The last exponential, which is 
of the form exp[— i-fp -1 ] with 7 := 3ha 2 t/(32c 3 ), 
does not change the momentum of the wavefunction 
but affects the quasi-particle's position according to 
exp[— i7p _1 ]i exp[«7p _1 ] = x + jhp~ 2 . This expression 
is not very insightful, however, so that we exploit our as- 
sumption that the quasi-particle mode is initially local- 
ized around some momentum pq. We then can expand 

the operator p~ 1 as 1 /p - (p - p ) jp% + (p - p ) 2 /pi H 

so that we essentially get expressions which are linear or 
quadratic in p. The linear part shifts the position analo- 
gous to Eq. (|l^) . The term quadratic in p has the same 
structure as the operator p 2 /(2M) of the kinetic energy 
of an ordinary Schrodinger particle. Hence an acceler- 
ation introduces a momentum-dependent mass cx Po/a 2 
for sound waves in a BEC. 

Putting everything together by applying successively 
all exponentials of Eq. (^) to the quasi-particle mode we 
arrive at 
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This expression coincides with the Heisenberg-picture po- 
sition operator xh{—£) — exp[— id> + t]x exp[iLo + t] if the 
variable x on the r.h.s. of Eq. ( [l6| ) is replaced by the 
position operator x. 

In view of the discussions given above the result ( |l5| ) is 
easy to understand. The prefactor in curly brackets cor- 
responds to a position dependent rescaling of the wave- 
function, a combined effect of the first exponential in 
Eq. @ and the prefactors present in Eqs. ( pd| ) and (|l3|). 
The phase factor stems from the last exponential. In 
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Xii(—t) the first term is a result of the squeezing trans- 
formation (pi]). The second term describes an accelera- 
tion of the quasi-particle motion. Since the quasi-particle 
has to overcome the potential created by the accelerated 
ground-state, which according to Eq. (||) has the oppo- 
site sign as V itself, the acceleration is only half as large 
as for a usual particle (prefactor 1/4 instead of 1/2). The 
term proportional to ct corresponds to the sound prop- 
agation, where a correction oc a 2 jp\ to the velocity of 
sound appears because of the last exponential in Eq. (|^). 
The last term, which is nonlinear in a;, is caused by the 
rational transformation of Eq. (|l^) . 

It is also interesting to study the momentum opera- 
tor t) := exp[— io) + t]j5exp[ia) + t] in the Heisenberg 
picture, which is given by 



, at a 2 t 2 



a 2 t _ a 2 t o 
(17) 



The first term describes, for a < 0, a squeezing of the 
momentum. The second term has its origin in the ra- 
tional transformation (|l4|) which also produces a higher 
order correction to the first term. The last term is pro- 
duced by the fifth exponential in Eq. (||) . Because of this 
term pu{—t) is not a hermitean operator anymore. This 
is also a consequence of the asymmetric occurence of x 
and p in the definition of w+ . 

For an ordinary accelerated particle the momentum in 
the Heisenberg picture is given by p + Mat where the 
last term corresponds to the increase of the momentum 
due to the acceleration. Notably this term is absent in 
the case ( |l7| ) of a quasi-particle. Thus, the momentum 
of a quasi-particle is squeezed rather than accelerated by 
a weak external homogeneous force. 

This behaviour can be understood if we consider — lu + 
of Eq. (0) as a kind of Hamiltonian for a classical system, 
i.e., H — ~p[c 2 + ax] 1 / 2 (the minus sign arises because 
we are dealing with Xh(— t) instead of xu{i))- It then 
becomes obvious that the classical canonical momentum 
is not proportional to dx/dt anymore so that its time 
evolution is more complicated. This is reflected by the 
squeezing transformation. We also remark that the re- 
sults ( pf) and © also serve as a further justification 
for Eq. (|8|) since for Ti = they do agree with the so- 
lutions of the classical equations of motion, x = dH/dp 
and p = —dH/dx. 

Up to now we have analysed the time evolution of the 
left-moving wavepacket Ql{x,£) := (x\ exp[— iuj + t\ |Oq). 
The corresponding right-moving part Qn(x,t) :— 
(x\ exp[iili + t\ |0o) can be obtained by reversing the sign of 
t in the results (|l5|), (16), and Jl7|). It then becomes ob- 
vious that (for a > 0) the momentum of O/? is squeezed 
instead of its position (see Fig. 1). 

In conlusion we have demonstrated that a weak ho- 
mogeneous force acting on a BEC alters the motion of 



quasi-particle modes in a different way than for ordinary 
particles. The acceleration of the position is smaller since 
the mode has to overcome the additional potential pro- 
vided by the ground state. The momentum is squeezed 
rather than accelerated because the "Hamiltonian" ui + 
contains a term oc px which acts in the same way as a 
squeezing operator in quantum optics. 
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